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abstract 


The  relative  resolution  and  stability  of  conventional  and  new 
methods  of  spectral  estimation  are  investigated.  Techniques  considered  are 
conventional  beamsteer,  maximum  likelihood,  maximum  entropy  and  principle 
components.  These  techniques  are  evaluated  using  both  synthetic  data  and  real 
seismic  data  from  the  TFO  extended  short-period  array. 
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SECTION  I  . 

INTRODUCTION 

.  i 

The  goal  of  this  research  is  to  investigate  the  relative  resolution 
and  stability  of  conventional  and  new  methods  of  array  processing  or  spectral  ' 
estimation.  Resolution  as  intended  here  will  mean  the  directivity  of  the  process¬ 
ing  technique  or  the  ability  to  differentiate  two  signals  from  almost  the  samb  direc- 
tton.  Also,  we  will  investigate  resolution  from  the  point  of  view  of  the  technique's 
ability  to  detect  weak  signals, i.  e.  .  the  relative  noise  threshold  of  the  techniques, 
These  questions  will  be  studied  with  true  covariance  matrices  corresponding  to 

1  f“ed  “Umber  °f  Pla"e  WaVe  s*6nals  plus  noise.  Stability  is  investigated  by  gen¬ 
erating  random  column  vectors  VJf  Ul.  .  . .  N  from  , he  covariance  matrix  fl  and 
applying  the  techniques  to  the  conventional  estimate  n=  N-'  ?  V.  V.H  of  the  cov¬ 
ariance  matrix,  where  H  denotes  conjugate  transpose.  ,  *  ‘  , 

Two  Simplifications  of  the  general  problem  have  been  made  so  that 
the  computations  become  reasonable.  First,  w,  will  conduct  the  analysis  at  a 
given  frequency  so  that  the  matrix  dimension  is  C  by  C  where  C  is  the  number  of  . 
sensors.  Second,  computation  for  the  maximum  entropy  speqtrum  is  greatly 

simplified  by  using  a  line  array  of  equally  spaced  sensors,  i.  e.  ,  we  avoid  iterative 
methods  for  spectral  estimation. 

I  ' 

The  techniques  to  be  considered  are  conventional  beamsteer  BS 
maximum-likelihood  (unbiased  minimum  mean  square  error)  ML,  maximum 
entropy  (K-Line  or  Burg  technique),  and  principal  components  or  eigenvalue  ’  '' 

analysis.  Exact  mathematical  description  of  these  techniques  will  be  given  in  the 
following  section.  1 

The  method  of  principal  components  differs  somewhat  in  philosophv 
from  the  other  techniques  and  i,  is  not  directly  a  method  for  spectral  estimation  ' 


1-1 
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t 

!  I 

or  array  processing.  ,  Justification  for  its  inclusion  lies  in  its  complete  invar¬ 
iance  to  random  noise  and  the  tendency  to  isolate  plane  wave  signals  on  eigen¬ 
vectors  ordered  by  signal  strength.  The  procedure  is  to  compute  the  eigenvectors 
and  then  compute  the  eigenvector  response  function  (or  eigenspectra),  which  is 
the  absolute  square  value  of  the  dot  product  of  the  eigenvector  and  a  "look  direc¬ 
tion"  vector.  The  direction  at  the  maximum  of  the  successive  eigenspectra  will 
be  shown  empirically  to  correspond  to  the  directions  of  the  plane  wave  signals 
specified  in  the  covariance  matrix.  Mathematically  QZ  =  Az  must  be  satisfied, 
where  Z,  and  A  are  an  eigenvector  and  eigenvalue  off}.  Thus  (f}+01)  Z=  (A+/3)Z 
for  any  constant  $,  so  that  the  eigenvectors  of  f}  are  also  the  eigenvectors  of 
n+01  and  the  technique  is  therefore^  expected  to  be  insensitive  to  random  noise. 

There  has  been  considerable  previous  theoretical  treatment  of  the 
principal  components  method  but  very  little  actual  application  of  the  technique  to 
array  processing.  The  classical  multivariate  statistics  approach  is  treated  by 
Anderson  (1958)  and  the  interpretation  of  these  results  in  terms  of  multiple  station¬ 
ary , time  series  is  discussed  by  N.  R.  Goodman  (1967).  The  actual  stimulation  for 
our  consideration  of  principal  components  in  the  present  effort  came  from  a  paper 
by  N.  L.  Owsley  (1971)  where  he  suggested  that  plane  wave  signals  tend  to  be  iso¬ 
lated  on  individual  eigenvectors  and  ordered  according  to  strength. 


i 

i 


i 
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SECTION  II 


/ 

MATHEMATICAL  DESCRIPTION  OF  TECHNIQUES  AND 
RESOLUTION  RESULTS 


A.  GENERAL 

A  synthetic  or  theoretical  cross-power  matrix  was  computed  and 
used  to  evaluate  the  following  four  types  of  spectra  estimates: 

•  Beamsteer  (BS) 

•  Maximum- Likelihood  (ML) 

•  Markoff  or  K-Line  (MK) 

•  Eigenvalue -Vector  (EG) 


The  cross-power  matrices,  fl  were  generated  from  special  signal 
vectors,  s^,  where  j  =  1,  2,  ....  N,  according  to  the  equation: 

N 

L,  a.s.s.  +  j8l 
jl;l  J  J  J 

where  the  a  's  and  $  are  the  relative  weights  for  the  individual  signals  and  ran¬ 
dom  noise  component  respectively,  and  H  consists  of  the  operation  of  transposi¬ 
tion  and  conjugation. 


The  elements  of  the  signal  vector  are  defined  as: 


s.  (/)  = 

J 


e-i27TfAt.  (£) 


,  1=  1,  2, 


C 


where  C  is  the  number  of  elements  or  array  sensors  to  be  used  in  the  evaluation, 
f  is  frequency,  and  At.  (/)  is  the  time  delay  for  the  /h  channel  in  the  jth  look 
direction.  The  array  can  be  three-dimensional,  (X,  Y,Z).  The  signal  vectors  are 
specified  by  an  azimuth  angle  (clockwise  direction  from  positive  Y-ordinate  in 
X-Y  plane),  elevation  angle  (angle  of  incidence),  and  a  propagation  velocity. 
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Power,  P,  is  computed  for  a  matrix  generated  at  a  selected  fre¬ 
quency  as  a  function  of  the  "look  direction"  (that  is,  the  direction  having  an 
azimuth  angle  of  0  )  for  each  of  the  previously  listed  techniques.  The  "look  dir¬ 
ection"  vector,  S($ )  is  defined  similarly  to  the  signal  vector  except  the  angle  of 
elevation  is  assumed  to  be  zero.  Starting  with  0=0  and  increasing  0  in  even 
increments  until  0=  27T  provides  a  complete  azimuth-power  spectrum. 

B.  MATHEMATICAL  DESCRIPTION  OF  THE  TECHNIQUES 

The  following  is  a  brief  description  of  the  procedure  used  to  cal¬ 
culate  these  spectra. 


1. 


Beamsteer  Spectra, 


PBS<«> 


The  beamsteer  spectra  are  the  response  of  the  beamed  (in  the  "look 
direction",  0  )  array  to  the  data  matrix  Q.  Effectively  the  beamed  array  has  the 
response  of  a  filter  vector  b  defined  as: 


b  =  S  (0)/  C 

so  that 

PBS(0)  =  bHfJb 

By  the  definition  of  fl  the  spectra  are  actually  computed  by  the  equation: 


PBS^ 


N 


+  fi/C 


Maximum-likelihood,  PML(0) 

The  maximum-likelihood  filter  M(0),  is  defined  by 
m(0)  =  n'1  S/(SH  A"1  S) 
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Therefore  the  maximum -likelihood  spectrum  is  given  by: 


PML<«>  =  l/js^JV1  S  (0)] 


3.  Maximum  Entropy,  Px,r,(0) 

MK 

The  technique  as  programmed  is  applicable  only  to  an  array  con¬ 
sisting  of  elements  with  equal  spacing  Ax  along  a  single  line  or  several  parallel 
lines.  In  the  case  of  several  lines  of  sensors  the  cross-power  matrix  for  each 
line  is  stacked  to  form  a  single  matrix  dimensioned  by  the  size  of  the  largest 
individual  line  matrix. 


The  inverse  of  is  computed  and  prediction  error  filters, 

are  determined.  The  azimuthal  power  spectra  is  then  computed  according  to  the 
following: 


=  2  p  Ax 


C 

L 

£=i 


-i27 T  kA  x 


-2 


where  p  is  the  prediction  error  and  k  is  the  wave  number,  which  is  given  by: 
k  =  f  cos  0  /V. 

where  V  is  the  velocity  of  propagation. 

4.  Eigenspectra,  P  E(-.{0) 

The  matrix,  fl  ,  can  be  represented  as 


n=  E 

i=i 


H 


7  h 8 i 


where  A^  and  g^  are  the  eigenvalues  and  eigenvectors  respectively  of  the  matrix 
0.  Having  developed  the  matrix,  fi  ,  from  the  specified  signal  vectors,  the 
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eigenvalues  and  eigenvectors  are  how  extracted  from  O  and  compared  to  the 
original  weighting  values  a.  and  signal  vectors  s  The  response  of  the  eigen¬ 
vectors  associated  with  the  eigenvalues  are  computed  as  the  eigenspectra,  P 
where  EG/ 


*  i=  1. 


C. 


The  ^  '8  are  arranged  in  decreasing  order  of  magnitude. 

C.  COMPARATIVE  RESOLUTION  CAPABILITIES  FOR  THE  SYNTHETIC 
DAT  A 

The  resolution  and  accuracy  of  each  of  the  spectral  estimates  were 
evaluated  from  a  cross-power  matrix  at  3.  0  Hz.  developed  from  seven  signal 
vectors  and  1  per  cent  random  noise.  The  specifications  for  the  seven  signals 
are  presented  in  Table  H-l;  note  that  the  signals  are  progressively  smaller  and 
that  Signal  7  is  50  dB  down  from  Signal  1. 


TABLE  II- 1 

SIGNAL  MODELS  FOR  SPECTRAL  EVALUATION 


Signal 

Number 

Velocity 

Azimuth 

Elevation 

Weighting 

Relative 
Power  (dB) 

1 

10  km/ sec 

0.  0 

0.  0 

10,000 

0.  0 

2 

!! 

40.  0 

0.  0 

1,000 

o 

o’ 

1 

3 

ft 

320.  0 

0.  0 

.100 

-20.  0 

4 

! ! 

340.  0 

0.  0 

.  0100 

-30.  0 

5 

1  1 

300.  0 

0.  0 

.  0050 

-33.  0 

6 

If 

20.  0 

0.  0 

.  0010 

o 

o’ 

■ 

7 

11 

60.  0 

0.  0 

.  0001 

-50.  0 
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The  array  consisted  of  2^0  equally  spaced  elements  at  1  km  inter- 

:  '  I 

vals  along  the  X-ordinatei  Thus  a  value  of  9=  0  corresponds  to  a  look  direction 
perpendicular  to  the  array.  ■  ! 

i 

!  I 

The  spectral  estimates  obtained  by  the  BS,  ML  and  MK  methods 

i  1  1 

are  shown  in  Figure  II- 1.  The  exact  location  of  the  signals  are  marked  by  vertical 
lines  at  the  top  of  the  figure.  The  MK  spectrum  is  plotted  by  scaling  its  maximum 
value  to  that  of  the  ML  spectrum.  Azimuth  is  plotted  over  360°  for  the  line  array  ’ 
so  the  interval  90  to  270  is  the  alias  of  the  remainder  of  th,e  plot.  An  analysis 
of  the  figure  shows  that  only  the  first  two  signals  (1  and  2),are  resolved  by  the  BS 
spectra.  Signals  1-5  were  resolved  by  both  the  ML  and  MK  spectra,  with  the  MK 

1  *  '  *  '  j 

spectra  exhibiting  much  sharper  and  stronger  peaks  than  the  ML  spectra. 

None  of  the  three  methods  were  able  to  resolve  the  weak  signal  6 
which  was  positioned  between  the  two  strpng  signals  1  and  2.  However,  an  even 
weaker  signal  (No.  7)  seems  to  have  been  resolved  by  the  MK  estimate.  To  det¬ 
ermine  if  signal  7  was  being  resolved  or  if  the  peak  is  just  a  sidelobe  effect,  a- 
r.other  spectral  computation  run  was  made  in  which  signal  7  was  moved  further 
away  from  signal  2.  As  shown  in  Figure  II-2  the  peak  did  not  move  Ivith  the  signal 
location  indicating  sidelobe  effects  and  the  inability  of  all  th^ee  techniques  to  re¬ 
solve  this  signal  energy.  •  1  ! 

s  '  '  i '  1 

In  the  second  run,  the  weighting  of  signal  6  was  increased  by  7  dB, 
leaving  its  level  33  dB  down  from  the  adjacent  signal  1.  In  this  case  (Figure  H-2) 
both  the  ML  and  MK  spectra  evidently  resolved  the  signal  pnergy.  ;  At  the  same 
time,  signal  5  was  decreased  by  7  dB  to  40  dB  down  from  signal  1.  Here  the  MK 
spectra  still  resolved  the  energy  but  the  Ml  spectra  tends  to  lose  the  energy  peak. 

i 

EG  spectra  Were  computed  from  the  cross-power  matrices  on  which 
the  previous  three  techniques  were  evaluated.  For  the  first  case,  seven  eigen¬ 
vectors  were  extracted  from  the  'matrix  and  spectral  estimates  made  from  each. 
Arranging  the  EG  spectra  in  order  based  on  the  decreasing  value  of  the  eigen¬ 
values,  a  one-to-one  correspondence  of  the  dominatihg  peak  in  each  of  the  spectra 

I 
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Igure  n-l «  spectral  Comparison  -  Seven  Signals 


to  the  signals  1  through  7  was  made, (Figure  II-3).  The  resolution  of  signals  6 
and  7,  which  were  previously  unresolved  by  the  other  techniques,  is  extremely 
encouraging.  Similarly,  the,  EG  spectra  from  the  second  run,  resolved  all  of  the 
signals.  Spectra  from  the  additional  eighth  and  ninth  eigenvectors  computed  in 
this  run  showed  many  lobes  within  6  dB  of  each  other  indicating  that  no  significant 
signal  energy  remained  in  the  matrix  after  the  extraction  of  the  first  seven  eigen¬ 
vectors.  (We  have  chosen  to  require  the  EG  spectra  to  have  a  peak  6  dB  above 
the  general  sidelobe  level  before  a  signal  component  is  "identified”.  ) 


To  further  test  the  EG  spectra  method,  a  set  of  signals  were 
chosen  in  which  several  pairs  of  the  signals  were  picked  to  be  highly  non-orthogonal 
(i.  e. ,  the  direction  of  the  second  signal  of  the  pair  was  in  the  direction  of  the 
maximum  sidelobe  of  the  first  signal). 


/.•> 


The  signal  models  in  this  case  are  listed  in  Table  11-2. 


TABLE  II-2 


4 


SIGNAL  MODELS  SELECTED  FOR  NON- OR THOGONAIiT Y 


Signal 

Number 

Velocity 

Azimuth 

Elevation 

Weighting 

Relative 
Power  (db) 

1 

10  km/  sec. 

o 

• 

o 

0.  0 

10. oOoo 

0.  0 

2 

II 

14.  0 

o 

• 

o 

1.  0000 

-10.  0 

3 

It 

29.0 

o 

• 

o 

.  1000 

-20.  0 

4 

II 

319.0 

0.  0 

.  0100 

-30.  0 

5 

II 

46.  0 

o 

• 

o 

.  0050 

-33.  0 

6 

11 

305.  0 

0.  0 

.  0010 

-40.  0 

7 

II 

73.  0 

~ I 

o 

• 

o 

.  0001 

-50.0 

t 


$ 

\ 
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Results  from  these  signal  models  were  similar  to  those  of  the 
previous  two  rune.  That  is,  the  BS  spectra  could  resolve  only  the  first  two  sig¬ 
nals,  ML  and  MK  spectra  resolved  the  first  5  signals,  and  the  EG  spectra  re¬ 
solved  all  seveu  signals. 

The  effects  of  isotropic  noise  on  the  ML,  MK  and  EG  spectral 
estimates  was  evaluated  using  the  signals  presented  in  Table  II-2,  (Figures  II-4 
and  H-5).  The  previously  discussed  run  did  not  contain  any  isotropic  noise.  The 
addition  of  1%  isotropic  noise  had  a  negligible  effect  on  the  EG  spectra.  The  most 
significant  effect  occurred  in  the  ML  spectra  where  the  "whole"  noise  level  in¬ 
creased  3  to  5  dB,  almost  obscuring  signal  5. 

The  addition  of  10  percent  isotropic  noise  prevented  both  the  ML 
and  MK  spectral  estimates  from  resolving  signals  4  and  5.  However,  the  EG 
spectra  still  resolved  all  seven  of  the  signals,  although  the  azimuth  of  the  weakest 
signal  was  biased  (Figure  II-5).  The  eigenspectra  for  the  stronger  signals  (1-5) 
were  not  significantly  effected  by  the  isotropic  noise.  The  sidelobes  of  the  eigen¬ 
spectra  corresponding  to  the  weaker  signals  increased  approximately  5  dB  due  to 
the  isotropic  noise,  however  the  signal  peaks  were  still  predominant  enough  to 
allow  detection  of  signals. 

Two  final  data  matrices  were  constructed  using  the  parameters 
m  Table  II-3  in  order  to  investigate  signals  extremely  close  in  azimuth  and  signals 
of  near  equal  strength.  The  signal  vectors  were  not  normalized  in  the  program 
so  the  approximate  eigenvalues  (assuming  independence  of  the  input  signal  vectors) 
are  20  times  the  input  signal  strength.,  Calculated  eigenvalues  are  shown  in  the 
table  and  are  similar  to  the  approximate  values,  given  in  parentheses. 

Signals  near  in  azimuth  were  grouped  into  azimuthal  pairs  separated 
by  only  a  few  degrees  and  were  selected  so  that  all  the  pairs  were  in  the  main  lobe 
of  the  beamsteer  response  pattern.  Spectra  obtained  from  the  BS,  ML  and  MK 
techniques  are  shown  in  Figure  II-6a,  for  the  EG  technique  in  Figure  II-6b.  The 
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Figure  II-4.  Effect  of  Isotropic  Noise  on  Spectra 
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Figure  II-6b.  Eigenspectra  For  Signals  Near  In  Azimuth 


TABLE  II- 3 

SIGNAL  OF  NEAR  STRENGTH  AND  AZIMUTH 


Near  Azimuth 

Near  Strength 

Signal 

Number 

Azimuth 

Strength 

Eigenvalue 

Azimuth 

Strength 

Eigenvalue 

1 

176 

10 

(200) 

0 

10 

(200)  209.6 

2 

185 

5 

(100) 

14 

5 

(100)  98.8 

3 

232 

1 

(20) 

29 

3 

(60)  56.3 

4 

226 

.  5 

(10) 

319 

2 

(40)  40. 2 

5 

130 

.  1 

(2) 

46 

1.5 

(30)  30.  0 

6 

132 

.  05 

(1) 

305 

1.2 

(24)  23.7 

7 

73 

1.  1 

(22)  19.0 

MK  method  has  the  best  capability  because  it  resolves  signals  3  and  4  as  well  as 
signals  1  and  2.  The  ML  and  EG  techniques  resolve  signals  1  and  2  only,  the  BS 
technique  does  not  resolve  any  pair. 

Results  for  the  near  equal  strength  data  were  the  same  as  those 
obtained  previously  and  the  figures  are  not  shown. 
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SECTION  III 
STABILITY  RESULTS 

The  spectral  estimation  techniques  evaluated  for  resolution  in  the 
preceeding  section  on  the  "true"  matrix  for  the  set  of  specified  signal  vectors  as 
defined  by: 

N 

n=  X  a.  s,  sH  +  51 
j=l  J  J  ' 

were  also  evaluated  for  stability  on  an  "estimate",  ,  of  the  true  matrix  f)  de¬ 
fined  as 


n.+J  =  yn.  +  a  -y)  x.H 

where  y  is  an  exponential  decay  factor  and  X ^  is  the  ith  generated  transform 
vector  based  on  the  set  of  specified  signals.  That  is, 

N  1/7 

X.  =  E  «.  s.  +  B  72  V. 

1  j=l  J  J 

where  c  is  an  independent  scalar  gaussian  deviate  such  that  E  U  «  =  a. 

j  L  j  j  j  j 

and  TJ.  is  a  normally  distributed  random  vector  with  zero  mean  and  identity  cov¬ 
ariance. 

An  estimate  of  the  matrix  could  be  obtained  by  an  average  of  L 
matrices,  each  generated  from  a  transform  vector,  x.  ,  as  previously  defined. 
That  is, 

ai  ■  | 
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The  stability  of  the  spectral  estimation  techniques  will  be  a  function  of  the  decay 
factor,  y .  The  number  M  {equivalent  to  L  for  successive  updates  of  Q.)  is  re¬ 
lated  to  y  by  the  equation: 

M  =  (i  +y)  /  (i  -y) 

The  first  stability  test  run  was  conducted  on  an  "estimated"  matrix 
computed  as  outlined  above  from  the  set  of  three  signals  described  in  Table  111-1, 
w^t^1  y=.95  and  the  number  of  updates  between  spectral  estimates  M  =  40. 


TABLE  III- 1 

FIRST  SIGNAL  MODEL  FOR  STABILITY  STUDIES 


Signal 

Number 

Velocity 

Azimuth 

Elevation 

Weighting 

1 

10  km/ sec 

240.  0 

0.  0 

2 

2 

II 

92.0 

0.  0 

1 

3 

II 

150.  0 

0.0 

.  5 

Starting  with  an  initial  matrix  generated  from  the  average  of  20 
aata  vectorSi  a  set  of  spectral  estimations  was  made  from  the  "estimated"  matrix 
after  40,  80,  and  120  updates.  The  spectra  after  40,  80,  and  120  updates  are 
shown  in  Figures  III- la  through  III-lc  with  the  BS,  ML  (which  is  plotted  28  dB 
down)  and  MK  spectra  at  the  top  and  the  three  EG  spectra  at  the  bottom.  The  MK 
spectra  were  developed  from  5  point  prediction  filters  thereby  resulting  in  some¬ 
what  lower  resolution  than  observed  in  the  previous  study.  The  stability  of  all 
estimation  techniques  is  extremely  good  for  this  case,  the  three  estimates  {for  a 
given  technique)  are  virtually  identical. 
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Figure  III- la.  Stability  of  Spectral  Estimate 
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Figure  III-lc.  Stability  of  Spectral  Estimates,  M=40 
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A  second  case  was  run  to  further  evaluate  the  stability  of  the 
estimation  in  which  the  decay  rate  was  increased  by  lowering  y  to  .  9  .  Seven 
signal  vectors  as  described  in  Table  III-2  were  used  in  this  case.  The  number 
of  updates,  M,  between  successive  spectra  computations  was  19. 

TABLE  III -2 

SECOND  SIGNAL  MODEL  FOR  STABILITY  STUDIES 


Signal 

Number 

Velocity 

Az  imuth 

Elevation 

Weighting 

1 

1 0  km/ sec 

180.  0 

0.  0 

10.  000 

2 

f  1 

240.0 

0.  0 

2.  000 

3 

ff 

92.0 

0.  0 

1.  000 

4 

If 

150.  0 

0.  0 

.500 

5 

120.0 

0.  0 

.  100 

6 

If 

210.  0 

0.  0 

.  010 

7 

ff 

270.  0 

0.  0 

.  001 

The  ML  and  MK  spectra  computed  after  19,  38,  and  57  updates 
are  presented  in  Figure  III-2.  For  the  MK  spectra,  a  10  point  prediction  error 
filter  was  used.  The  stability  of  the  estimates  for  the  signals  1  to  4  appear  to  be 
very  good  for  all  processes  with  one  exception;  the  MK  spectra  was  slightly  un¬ 
stable  in  the  ar°a  of  signal  3. 

Signal  5  was  resolved  only  by  the  ML  estimate  and  the  stability 
seems  to  be  fair.  Signal  6  was  resolved  by  both  the  ML  and  MK  estimates,  how¬ 
ever  the  MK  line  is  rather  unstable.  The  seventh  signal  was  not  resolved  by  any 
of  the  three  techniques. 


Ill- 6 


Similarly,  EG  spectra  were  computed  for  the  "estimated"  matrices 
(at  the  end  of  19,  38.  and  57  updates).  The  spectra  computed  from  the  first  eight 
extracted  eigenvectors  of  each  matrix  are  shown  in  Figure  111-3.  The  stability 
of  the  spectra  in  the  area  of  the  signal  appears  to  be  very  good. 


The  first  six  signals  are  detectable  in  all  sets  of  the  EG  spectra. 

Note  that  signals  2,  3,  and  4,  whose  specified  power  levels  relative  to  the  strong¬ 
est  signal  are  -7,  -10.  and  -13  dB  respectively,  appear  to  be  somewhat  mixed  in 
the  second,  third,  and  fourth  EG  spectra.  That  is.  spectrum  No.  2  contains  pri¬ 
marily  signal  2  with  an  indication  of  signal  3  and  spectrum  No.  3  contains  primarily 
signal  4  with  an  indication  of  both  signals  2  and  3.  However,  because  these  "over¬ 
lapped"  signals  are  6  dB  or  more  below  the  "target"  signal,  it  is  still  possible  to 

unambiguously  identify  signal  S.  from  EG  spectrum  No.  S 

1  ’  i' 

Although  spectrum  No.  7  consistently  peaked  up  at  the  location  of 
signal  7,  the  signal  would  not  have  been  detectable  from  the  computed  spectrum, 
because  the  sidelobes  are  as  large  as  the  actual  peak.  However,  if  spectrum  No.  7 
were  averaged  over  time,  that  is  after  updates  19,  38,  and  57.  etc. .  signal  7  might 

be  detectable  because  of  the  unstable  nature  of  the  spectra  outside  the  area  of  signal 
7  relative  to  the  stability  in  that  area. 
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SECTION  IV 

TFO  EIGENSPECTRA  ANALYSIS 


Seismic  noise  data  from  the  extended  short  period  array  located 
at  TFO  was  analyzed  under  previous  contract  F33657-70-C-0I  00  using  a  high 
resolution  (HR)  spectral  technique.  This  HR  technique  is  of  the  form 
l/s<#)  fl-  S(«)  and  is  described  by  Texas  Instruments  Incorporated  (1970).  The 
technique  generally  gives  reasonable  wavenumber  spectra,  but  occasionally  has 
spurious  secondary  peaks.  The  HR  estimates  shown  here  are  thought  to  be  satis- 
factory;  they  were  chosen  for  comparison  because  they  were  readily  available. 
One  of  the  noise  samples  used  in  the  analysis  (designated  as  "winter"  noise)  was 
selected  for  the  evaluation  of  EG  spectral  analysis. 


A.  REVIEW  OF  PREVIOUS  ANALYSIS 

A  significant  increase  in  noise  power  (approximately  12  dB)  was 
observed  for  the  "winter"  sample  and  was  due  partially  to  intense  storm  activity 
off  the  west  coast  of  Alaska.  At  low  frequencies  (.  12  Hz)  a  significant  amount  of 
surface  mode  noise  generated  by  the  Alaskan  storm  is  evident  (Figure  IV- 1).  The 
surface  mode  contribution  from  this  storm  decreases  at  higher  frequencies  and 
eventually  a  source  to  the  Northeast  becomes  the  predominant  surface  mode  noise, 
(at  about  .  70  Hz).  At  this  and  higher  frequencies,  the  P-wave  noise  dominates 

the  spectra.  Analysis  of  the  P-wave  noise  at  frequencies  around  .  5  Hz  revealed 
several  sources  (Figure  IV-2). 


EIGENSPECTRA  COMPUTATIONS 


Based  on  the  previous  noise  analysis  of  the  "winter"  noise  sample, 
the  crosspower  spectrum  matrix  computed  at  the  frequency  of  .  284  Hr  was  seleci 
ed  for  the  eigenspectra  analysis.  At  this  frequency,  the  multi-directional  nature 
of  the  P-wave  energy  is  indicated  but  not  resolved.  Also  a  strong  surface  mode 
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noise  source  from  the  Northeast  and  a  weak  one  from  the  West  is  evident.  The 
HR  spectrum  computed  at  this  frequency  is  shown  in  Figure  IV-3;  note  the  iso¬ 
tropic  nature  of  all  surface  mode  noise. 

Ten  eigenvector -values  were  extracted  from  the  crosspower 
mat:  ix.  A  maximum  of  100  iterations  were  allowed  in  the  algorithm  to  extract 
the  eigenvalues.  If  the  criterion  was  not  satisfied  at  that  time  the  vector  obtain¬ 
ed  to  that  point  was  accepted  and  extracted  from  the  matrix.  Table  IV-1  pre¬ 
sents  the  number  of  iterations  required  and  resulting  eigenvalue  for  each  of  the 
eigenvectors. 


TABLE  IV-1 

NUMBER  OF  ITERATIONS  AND  RESULTING  EIGENVALUE 


Eigenvector 

Number 

Number  of 
Iterations 

Eigenvalue 

1 

19 

11. 322 

2 

23 

6.  005 

3 

22 

3.654 

4 

33 

2.  065 

5 

64 

1.488 

6 

66 

1.240 

7 

63 

1.  061 

8 

100 

.  902 

9 

66 

.828 

10 

100 

.699 

C.  WAVENUMBER  ANALYSIS  OF  EIGENVECTORS 

The  wavenumber  spectrum  of  each  eigenvector  was  computed  and 
results  are  compared  with  the  HR  spectrum  in  a  previous  analysis.  The  spectrum 
from  the  first  eigenvector  is  shown  in  Figure  IV-4,  The  spectrum  shows  a  very 
strong  peak  which  corresponds  to  the  dominating  P-wave  energy  from  the  Alaska 
storm  direction. 
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Figure  IV-4.  TFO  Eigenspectrum  1 
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The  spectra  from  the  next  two  eigenvectors  (Figure  IV-5)  show 
two  other  distinct  P-wave  contributions  to  the  noise  field  which  closely  agree 
with  the  P-wave  sources  observed  at  higher  frequencies  in  the  previous  analysis 
(See  Figure  IV-2), 

The  spectra  from  eigenvectors  No.  4  and  No.  5  (Figure  IV- 6) 
show  peaks  corresponding  to  surface  mode  energy,  propagating  from  the  North¬ 
east  of  the  array.  The  peak  powers  in  the  two  spectra  differ  by  about  10  degrees, 
indicating  the  presence  of  surface  mode  sources  to  the  Northeast  (or  possibly 
multipath  propagation  from  a  single  source).  Previous  HR  spectra  at  the  fre¬ 
quency  of  0.  33  Hz  (Figure  IV-7)  also  indicated  two  surface  mode  peaks  in  this 
direction.  This  noise  may  be  generated  by  low  pressure  activity  in  the  area  of 
Greenland  and  Eastern  coast  of  the  United  States. 

The  next  three  EG  spectra  (Figure  IV-8)  show  peaks  corresponding 
to  surface  mode  noise  from  the  west;  these  also  are  suggested  in  the  HR  spectrum 
of  Figure  IV-7.  Note  that  these  components  are  not  certain  because  the  EG  spectra 
are  showing  significant  sidelobes. 

The  spectrum  for  the  last  eigenvector  extracted  (Figure  IV-9)  has 
many  peaks  in  the  surface  mode  region,  suggesting  that  the  remaining  propagat¬ 
ing  noise  is  isotropic. 
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Figure  IV-5.  TFO  Eigenspectra  2,  3 
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Figure  IV-7.  TFO  Wavenumber  Spectrum1  at  0.  33  Hz 
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Figure  IV-8.  TFO  Eigenspectra  6,  7,  8 
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Figure  IV-9.  TFO  Eigenspectrum  9 
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SECTION  V 
CONCLUSIONS 


Four  wavenumber  spectral  estimation  techniques,  beamsteer  (BS), 
maximum  likelihood  (ML),  maximum  entropy  (MK)  and  eigenspectral  analysis  (EG) 
have  been  compared  in  terms  o£  resolution  and  stability  using  synthetic  covariance 
matrices  generated  for  a  linear  array.  In  addition,  the  EG  technique,  which  has 
not  been  applied  previously  in  seismology  was  used  or.  real  seismic  data  and  com- 

pared  with  the  ML  technique  applied  to  the  same  data.  Major  conclusions  are 
given  below: 


1.  For  the  synthetic  data  the  EG  technique  recovered  weak  signals 
best;  in  certain  cases  it  was  capable  of  recovering  signals  50  dB  below  the  largest 
signal  in  the  covariance  matrix.  The  ML  and  MK  techniques  could  recover  signals 
30  dB  below  the  largest  signal.  The  better  performance  by  the  EG  technique  re- 
suited  because  it  is  insensitive  to  random  noise. 

2.  The  MK  technique  showed  the  best  capability  to  resolve  two  signals 
having  similar  azimuths.  For  the  synthetic  example  used,  the  MK  technique  could 
resolve  two  signals  6°apart,  but  could  not  resolve  two  signals  2°  apart.  The  ML 

and  EG  techniques  could  resolve  signals  9°  but  not  6°  apart;  the  BS  technique  could 
not  resolve  signals  9°  apart. 


3.  The  stability  of  all  estimates  appeared  to  be  good.  The  BS  was  the 
most  stable,  the  EG  (at  least  in  the  vicinity  of  the  signal  directions)  seemed  to 
have  the  next  best  stability,  and  the  ML  and  MK  techniques  appeared  to  be  similar. 


4.  Application  of  the  EG  technique  to  real  seismic  data  gave  interpret¬ 
able  results  consistent  with  HR  analysis  of  the  same  data.  In  fact,  the  EG  tech¬ 
nique  seemed  to  give  better  resolution  of  the  propagating  noise  components. 
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5.  The  EG  technique  appears  to  be  a  very  valuable  tool  for  array 
analysis!  its  insensitivity  to  random  noise  is  a  very  important  property  because 
real  seismic  data  usuaUy  has  a  significant  random  noise  component  which  can 
limit  the  ML  and  MK  techniques.  Possible  applications  of  this  technique  include 
separation  of  overlapping  events  and  array  detection. 
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